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Abstract 

A high precision study of the process 77 ^ pp has been performed using a data 
sample of 89 fb“^ collected with the Belle detector at the KEKB e"’'e“ collider. 
The cross section of pp production has been measured at two-photon center-of- 
mass (c.m.) energies between 2.025 and 4.0 GeV and in the c.m. angular range of 
I cos 0*1 < 0.6. Production of 77 ^ 7 c ^ PP is observed and the product of the 
two-photon width of the pc and its branching ratio to pp is determined. 

PACS numbers: 12.38.Qk, 13.60.Rj, 14.40.Gx 
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1 Introduction 


Two-photon collisions provide a clean environment for baryon pair produc¬ 
tion and such events can be produced in great abundance at a high luminosity 
electron-positron collider. Accurate measurements of such processes, in par¬ 
ticular 77 — pp, is important to test existing theoretical predictions. 

General theories of hard exclusive processes in QCD [1,2] (see also [3] for 
a review) predict that the differential cross section for 77 — hih 2 at large 
energies and hxed c.m. angle {6*) has the form 

da/dt oc {as s ^ 00 ). ( 1 ) 


Here Uc is the number of elementary constituents participating in the hard 
interaction, f{6*) is a specihc function expressed via dehnite integrals over 
hadron wave functions, s is the square of the c.m. energy of the two-photon 
system, and t is the square of the four-momentum transfer from a photon to 
hadron. The first estimate of the cross section for 77 — pp was obtained in 
the three-quark picture {uc = 8 ) [4,5], using the proton wave function based 
on QCD sum rules [ 6 ]. Previous measurements [7,8,9,10,11] in the W^^{= ^/s) 
range between 2.5 and 3.0 GeV gave cross sections one order of magnitude 
larger than this expectation. To explain these experimental observations, var¬ 
ious model-dependent approaches were suggested. For example, in the diquark 
model [12,13,14] the proton is considered to be a quark-diquark system. In this 
case Uc = 6 and a diquark form factor is introduced, so that Eq. (1) becomes 
da/dt oc s~*\F\‘^, where F may depend on s. Asymptotically, F — f{6*)/s 
[12], and the behavior da/dt oc is recovered. These results exhibit bet¬ 
ter agreement with measurements of the absolute size of the cross section for 
= 2.5 - 3.0 GeV. 

Other approaches have been developed recently. The handbag model [15] has 
been developed for large momentum transfers, and the calculations have been 
applied at medium energies (IT77 > 2.55 GeV) with large uncertainty bands. In 
Ref. [16], the Veneziano model is applied in an unmodified form to the process, 
and fair agreement with data is obtained without adjustable parameters. 

Recently, the measured energy range for 77 —> pp has been extended to 
IF 77 = 4 GeV and above [10,11], but with very limited statistics for > 3 
GeV. Furthermore, pp 77 experiments give the cross section for the in¬ 
verse process at = 3.2 - 3.7 GeV [17]. To test QCD predictions, it is very 
important to improve the statistics at higher energies. Moreover, an accurate 
measurement with higher statistics for 77 —> pp is crucial in the study of 
the interactions involved. This paper presents the Belle measurement of the 
77 PP cross section for between 2.025 to 4.0 GeV and | cos 0*1 < 0 . 6 , 
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using a data sample corresponding to an integrated luminosity of 89 fb ^ 


2 Experimental apparatus and event selection 


Experimental data are recorded with the Belle detector [18] at KEKB [19], 
which is an asymmetric e“''e“ collider running at 10.58 GeV c.m. energy. In the 
laboratory frame, the direction of the positron beam is taken to dehne the — 
direction. For the analyses in this paper, the following Belle subsystems are of 
importance: the central drift chamber (GDC), the aerogel Gherenkov counters 
(AGG), the time-of-flight scintillation counters (TOE) and the Gsl electromag¬ 
netic calorimeter (EGL), all of which are located in a 1.5 T solenoidal magnetic 
held. The GDG measures the momenta of charged particles and provides par¬ 
ticle identihcation information by precise ( 6 %) dE/dx measurement, allowing 
separation of protons from other particles for momentum up to 1 GeV/c. The 
TOE measures the time of hight of particles with a 0.1 ns timing resolution, 
which is powerful for pjK separation for momentum up to 2 GeV/c. The EGL 
detects photons and is used to reject electrons by measuring the deposited en¬ 
ergy with a resolution of ge/E = 1.5% (2.0%) at 1 GeV (0.1 GeV). Using the 
number of photoelectrons observed, the AGG extends particle identihcation 
beyond that of the GDG dE/dx and TOE and is ehective in the suppression 
of highly relativistic and up to a momentum of 3.5 GeV/c. 

Through the process > e’''e“ 7 * 7 * —> e^e~pp, exclusive pp pairs are 

produced in quasi-real two-photon collisions, where the scattered e~^ and e~ 
are lost down the beam-pipe, and only the p and p can be detected. The 77 
axis thus approximates the beam direction in the e’''e“ c.m. frame. Gandidate 
events are searched for in a data stream where the sum of the magnitudes 
of momentum of all charged tracks and the total EGL energy are restricted 
to 6 GeV/c and 6 GeV, respectively. Events are required to have exactly two 
tracks of opposite charge satisfying the following conditions in the laboratory 
frame: pt > 0.35 GeV/c, dr <1 cm and \dz\ < 5 cm. Here pt is the transverse 
momentum and dr and dz are the radial and axial coordinates of the point of 
closest approach to the nominal collision point, respectively. Both tracks are 
required to hit the TOE counters. The invariant mass of the two tracks and 
the squared missing mass of the event, assuming the two tracks are massless, 
are required to be smaller than 4.5 GeV/c^ and larger than 2 GeV^/c^, re¬ 
spectively. A good transverse momentum balance in the e’''e“ c.m. frame is 
also required: |Sp)‘| = \p// + p//\ < 0.2 GeV/c, where p// and p)* denote the 
transverse momenta of the two tracks in that frame, with respect to the 
beam-axis. 

The selected events are dominated by 77 —> , n~^n~ and K^K~ up 

to this stage. Events with pp are separated from the others by a particle iden- 
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Fig. 1. Two-dimensional distribution of the normalized likelihood for pp identifi¬ 
cation, for the events passing all selection criteria except the cuts on the normal¬ 
ized likelihood indicated by the arrows. The dark parts show the events satisfying 
Ax = Tx /{Lp + Lk + + Lp + Le) > 0.01 for both tracks using the same x, where 
X is K, TT, p or e. 

tification (PID) algorithm, which is applied to each individual track under the 
following conditions: 1) the difference between the measured and the expected 
CDC dE/dx is less than 4 times the resolution: x^dE/dx = /dx)/adE/dxY < 

4^; 2) the ratio of the associated ECL energy to the momentum is less than 
0.9, which is only applied to the positively charged track; 3) the number of 
the photoelectrons in ACC counters associated with the track is less than 4, 
and this condition removes a large part of high-momentum and e^; 

4) the likelihoods for each particle assignment are combined to determine the 
normalized likelihood, \p = Lp/{Lp + Lk + Lt^ + Lp + L^), which has to be 
larger than 0.8 (Fig. 1). In these likelihoods, L = exp[—^{XdE/dx + Xr)] is 
calculated using information from the CDC {dE/dx) and the TOF (time of 
flight T). Here Xt = is the difference between the measured 

and the expected values for T, and aT is the timing resolution. A combined 
use of CDC and TOF allows p{p) separation from other particles, in particular 
K^, for momentum up to 2 GeV/c. 

For hFy.y = 2-4 GeV a total 36094 events survive all of the selection criteria. 
Their distribution in is shown in Fig. 2. A peak around 2.98 GeV can be 
identihed as the ?7c(2980) resonance [20]. A much narrower peak at 3.08 - 3.10 
GeV, corresponding to an excess of 26 ± 8 events relative to the continuum, 
could be attributed to backgrounds from radiative return to J/t/’, and the 
enhancement is in agreement with expectations based on this assumption [21]. 
Fitting the data from 2.6 to 3.7 GeV with a smooth (exponential of a hfth- 
order polynomial) function for the continuum, a Breit-Wigner function for 
the Pc and a Gaussian for the J/i/, a total pc yield of 156.9 ± 33.3 events 
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(GeV) (GeV) 

Fig. 2. distribution of events passing all the selection criteria: (a) = 2-4 

GeV; (b) a close-up view of (a) for = 2.82 - 3.26 GeV. 

is obtained. The statistical significance of the ? 7 c signal is 5.3 a, dehned as 
^—2 ln(Lo/T 5 ), where Ls and Lq denote the maximum likelihoods of the fits 
with and without a signal component, respectively. 


3 Monte Carlo simulation 


Monte Carlo samples for the following channels have been generated: e’''e“ —> 
e’''e“X where X is pp, K~^K~, e~^e~ and ppTi^. Hadron pair, 

lepton pair and 77 — ppn^ events are generated by the codes TREPS [22], 
AAFH [23] and GGLU [24], respectively. Event generation is followed by a 
detector simulation based on GEANT3 [25] and a trigger simulation. The 
selection criteria described in Section 2 are then applied to these Monte Carlo 
events. 

Because the acceptance depends on both and | cos 6* \, the signal (77 — 
pp) efficiencies are determined for a number of two-dimensional bins of the two 
variables. Other channels are generated for the study of residual backgrounds. 
Similar to 77 —> pp, the selection efficiencies for 77 — K^K~ and 77 — ppir^ 
are evaluated within each narrow bin. For the 77 — 7 r’'' 7 r“, and 

channels, realistic distributions are generated [23,26,27]. Events from those 
samples that survive the selection criteria are referred to as the expected 
residual backgrounds. 

From Monte Carlo simulation, the overall efficiency of 77 —pp for | cos0*| < 
0.1 ranges from ~ 3% at = 2 GeV to ~ 32% at = 4 GeV (Fig. 3). 
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Fig. 3. Overall detection efficiency of 77 pp as a function of and | cos 9* 

4 Measurement of the cross sections for 77 —pp 


Dividing | cos 6^*1 into bins of 0.1, W.y^ into bins of 25 MeV (for 2-3 GeV) 
and 100 MeV (for 3-4 GeV), the number of events selected from the data, 

I cos6^*1), is determined for each of the two-dimensional bins. The 
efficiency |cos0*|) is also evaluated from Monte Garlo simulation for 

each bin. The ratio of AiV to e is then converted to the differential cross 
section, according to the formula: 

da^^^ppiW^^) ^ AiV(l - f)/e 
dj cos9*\ Lint ^ AIV.,., A| cos0*| 


where / is the fraction of residual background in the data. Lint is the integrated 
luminosity and 7 ^?^ is the luminosity function. Here L^t = 88.96 fb“^, with 

a systematic uncertainty of 1.4%. The luminosity function 7 ^?^, as a function 
of IV,.,, is dehned by 




^e~''e pp 


nn I 




(IV 


dL^^iW. 


11) 


11) 


dW. 


dW. 


77 


77 


( 3 ) 


and is calculated by TREPS [22] using the equivalent photon approxima¬ 
tion method [21]. For the calculation of the luminosity function, the effects 
from longitudinal photons are neglected. For simulation in TREPS, the max¬ 
imum virtuality of each of the two photons, Q\ and Q 2 , is limited to 1 GeV^. 
Moreover, a form factor term is introduced for the high-Q^ suppression effect, 
{1 + {1 + . The systematic uncertainty of the luminosity 
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(GeV) 


Fig. 4. Measured cross sections for 77 ^ pp. For the Belle, CLEO [ 8 ] and VENUS 
[9] results, the error bars are purely statistical; while for OPAL [10] and L3 [11], 
both statistical and systematic uncertainties are included. Theoretical prediction 
curves shown are from [14] (diquark) and [4] (three-quark). 

function is estimated by comparing the kinematic distributions of the two- 
photon system for the events generated with TREPS to those from a QED 
calculation that includes all order diagrams [23]. Within the range = 2 
- 4 GeV agreement within 3 - 5% was reported [22,28]. 

The cross section a^^^ppiW^^) is obtained by a summation over ] cos 6 '*]: X] 
{da^^^ppiW^^)/d\ cos 6 '*]) A] cos 6 **], with the restriction of ] cos 6 **] < 0.6, due 
to polar angular coverage limits of the TOE system. The results are summa¬ 
rized in Table 1 and Fig. 4. The contribution from 77 —? 7 c —^ pp is included 
in these results. For the two cross section measurements in the lowest 
bins (2.025 - 2.075 GeV), the efficiencies are extremely small at larger ] cos 6 '*] 
and data are only available up to ] cos 6 '*] = 0.4 and 0.5, respectively. We thus 
fit a second-order polynomial function of cos^ 0 * to these differential cross sec¬ 
tions and arrive at a result by integrating the fit over ] cos0*] up to 0.6. For 
comparison, we also show in Fig. 4 the results from previous measurements 
[8,9,10,11]. 
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(GeV) 

cj(77 

pp) (nb) 

(GeV) 

CJ 

(77 

-^pp) 

(n 

b) 

2.025 - 2.050 

3.95 

± 

0.25 

± 

0.28 

2.650 - 2.675 

0.50 

± 

0.03 

± 

0.05 

2.050 - 2.075 

5.79 

± 

0.40 

± 

0.41 

2.675 - 2.700 

0.50 

± 

0.04 

± 

0.05 

2.075 - 2.100 

6.48 

± 

0.29 

± 

0.46 

2.700 - 2.725 

0.39 

± 

0.03 

± 

0.04 

2.100 - 2.125 

6.49 

± 

0.17 

± 

0.49 

2.725 - 2.750 

0.36 

± 

0.03 

± 

0.03 

2.125 - 2.150 

6.64 

± 

0.16 

± 

0.50 

2.750 - 2.775 

0.31 

± 

0.03 

± 

0.03 

2.150 - 2.175 

6.54 

± 

0.14 

± 

0.50 

2.775 - 2.800 

0.24 

± 

0.02 

± 

0.02 

2.175 - 2.200 

5.97 

± 

0.13 

± 

0.47 

2.800 - 2.825 

0.25 

± 

0.03 

± 

0.02 

2.200 - 2.225 

5.26 

± 

0.11 

± 

0.42 

2.825 - 2.850 

0.18 

± 

0.02 

± 

0.02 

2.225 - 2.250 

4.77 

± 

0.10 

± 

0.39 

2.850 - 2.875 

0.17 

± 

0.02 

± 

0.02 

2.250 - 2.275 

3.91 

± 

0.09 

± 

0.32 

2.875 - 2.900 

0.17 

± 

0.02 

± 

0.02 

2.275 - 2.300 

3.38 

± 

0.08 

± 

0.28 

2.900 - 2.925 

0.111 

± 

0.017 

± 

0.010 

2.300 - 2.325 

2.86 

± 

0.07 

± 

0.24 

2.925 - 2.950 

0.101 

± 

0.016 

± 

0.009 

2.325 - 2.350 

2.56 

± 

0.07 

± 

0.22 

2.950 - 2.975 

0.143 

± 

0.019 

± 

0.013 

2.350 - 2.375 

2.34 

± 

0.07 

± 

0.21 

2.975 - 3.000 

0.128 

± 

0.017 

± 

0.011 

2.375 - 2.400 

1.94 

± 

0.06 

± 

0.17 

3.000 - 3.100 

0.077 

± 

0.007 

± 

0.009 

2.400 - 2.425 

1.72 

± 

0.06 

± 

0.16 

3.100 - 3.200 

0.042 

± 

0.006 

± 

0.004 

2.425 - 2.450 

1.59 

± 

0.05 

± 

0.15 

3.200 - 3.300 

0.023 

± 

0.005 

± 

0.002 

2.450 - 2.475 

1.40 

± 

0.05 

± 

0.13 

3.300 - 3.400 

0.018 

± 

0.004 

± 

0.002 

2.475 - 2.500 

1.27 

± 

0.05 

± 

0.12 

3.400 - 3.500 

0.0092 ± 

0.0026 ± 

0.0009 

2.500 - 2.525 

1.21 

± 

0.05 

± 

0.11 

3.500 - 3.600 

0.0066 ± 

0.0021 ± 

0.0007 

2.525 - 2.550 

1.03 

± 

0.05 

± 

0.10 

3.600 - 3.700 

0.0090 ± 

0.0030 ± 

0.0010 

2.550 - 2.575 

0.88 

± 

0.05 

± 

0.08 

3.700 - 3.800 

0.0030 ± 

0.0019 ± 

0.0003 

2.575 - 2.600 

0.84 

± 

0.05 

± 

0.08 

3.800 - 3.900 

0.0045 ± 

0.0022 ± 

0.0005 

2.600 - 2.625 

0.74 

± 

0.04 

± 

0.07 

3.900 - 4.000 

0.0020 ± 

0.0012 ± 

0.0003 

2.625 - 2.650 

0.60 

± 

0.04 

± 

0.06 








Table 1 

Measured cross sections for 77 ^ pp (| cos0*| < 0.6). The first error is statistical 
and the second is systematic. 


Fig. 5 shows the angular dependence of the differential cross sections measured 
in 11 ranges of separately. The ? 7 c region (2.9 - 3.1 GeV) is skipped. For 
< 2.4 GeV, the differential cross sections decrease as | cos 0*1 increases; 
while for > 2.6 GeV, the opposite trend is observed. The transition 
occurs around = 2.5 GeV. Similar results are shown in Fig. 6 for three 
larger ranges of and are summarized in Table 2. For comparison, previous 
measurements [8,10,11] are also shown in Fig. 6. Further discussion is given 
in Section 6. 


Based on studies from Monte Garlo and data, residual backgrounds due to 
particle misidentihcation and non-exclusive events are subtracted from the 
data. Gorrections are based on AiV multiplied by (1 — /) as shown in Eq.(2). 
Gomplete details are given in Section 5, and the systematic errors are shown 
in Table 3. All measured cross sections and differential cross sections shown 
in this paper have been corrected in this way. The excess caused by the J/ijj 
background described at the end of Section 2 is estimated in each | cos 6* \ bin 
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yy 
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2.5 < W_< 2.6 GeV 
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2.7 < W_< 2.8 GeV 
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31 < W_< 3.5 GeV 3.5 < < 4.0 GeV 

yy yy 



Fig. 5. Measured differential cross sections in 11 ranges of as a function of 
I cos0*|. The error bars are statistical only. 


separately and subtracted from AiV before the other corrections above. The 
systematic uncertainty due to this subtraction is 8% for the measured cross 
section in the 3.0 - 3.1 GeV bin, taking into account the fluctuation of 
the estimated number of J/'ip. 
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2.0 < < 2.5 GeV 

7Y 



2.5 < W^<3.0 GeV 

yy 



3 < < 4 GeV 

yy 



Fig. 6. Measured differential cross sections as a function of | cos0*| for three ranges 
of The error bars are statistical only. 


COS0* 

d|cos6»*[ 

=2.075 - 2.5 GeV 

d[ cos 0*1 

=2.5 - 3.0 GeV 

=3 - 4 GeV 

0.0 - 0.1 

0.1 - 0.2 

0.2 - 0.3 

0.3 - 0.4 

0.4 - 0.5 

0.5 - 0.6 

10.28 ± 0.10 ± 0.29 

9.00 ± 0.10 ± 0.26 

7.20 ± 0.10 ± 0.27 

5.08 ± 0.10 ± 0.19 

3.79 ± 0.12 ± 0.17 
2.96 ± 0.17 ± 0.13 

0.56 ± 0.02 ± 0.02 

0.48 ± 0.02 ± 0.02 

0.53 ± 0.02 ± 0.04 

0.69 ± 0.03 ± 0.05 

0.94 ± 0.04 ± 0.03 
1.22 ± 0.05 ± 0.04 

0.010 ± 0.002 ± 0.002 

0.014 ± 0.002 ± 0.001 

0.013 ± 0.002 ± 0.001 

0.022 ± 0.003 ± 0.002 

0.038 ± 0.005 ± 0.003 

0.098 ± 0.010 ± 0.007 


Table 2 

Measured differential cross sections versus |cos0*| for different 1T^7 ranges. The 
first error is statistical and the second is systematic. 
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Source 

Systematic error (%) 

Integrated luminosity 

1.4 

Luminosity function 

3 - 5 

Trigger efficiency 

5 

Particle identification efficiency 

1 - 6 

Monte Carlo statistics 

1 - 3 

Particle misidentification background subtraction 

0 - 1 

Non-exclusive background subtraction 

2 - 12 

\'Lpl\ effect and residual non-exclusive backgrounds 

2 - 3 

J/t/j subtraction (IF 77 ~ - 3.1 GeV) 

8 

Possible backgrounds from radiative return 

1 

Total 

7- 14 


Table 3 

Systematic errors for the measured cross sections of 77 ^ pp. Some uncertainties 
are ITiy-y-dependent and shown as ranges. 

The total rjc yield, = 156.9 ± 33.3, can be converted to the product of 
the two-photon width of pc and the branching fraction of pc PP'- ^-yyiVc) x 
B{r]c pp) = e Lint = 7.20 ± 1.53(stat.)lg;?^(syst.) eV, 

using the luminosity function dL^^/dW^^ determined at the energy of the rjc 
mass {rririJ and the efficiency e from Monte Carlo. For the systematic error, 
effects from the uncertainties of the continuum background shape and the 
signal width are taken into account, in addition to all other sources listed in 
Table 3. The above result gives r.y.y(? 7 c) = 5.5±1.2(stat.)lo;g(syst.)±1.7(norm.) 
keV, where the last error comes from the branching fraction B{r]c —^ pp) 
uncertainty [20]. Since observations of the pc in the pp channel are scarce 
and suffer from low statistics, current measurements for B ( 77 ^ — pp) available 
in Ref. [20] are not very consistent with each other. Our result is the first 
measurement of T^^irjc) ^ B{rjf. ^ PP) in two-photon collisions and, together 
with the observation of the rjc in pp collisions in its 77 decay mode [29], will 
help to decrease the errors on both the pc two-photon width and branching 
fraction. 


5 Corrections and major sources of systematic error 


The accuracy of the Monte Carlo trigger efficiency has been checked from 
the two-track trigger, which requires at least two CDC tracks with an opening 
angle larger than 135°, two or more TOF hits as well as the ECL timing signal, 
using experimental events passing the high energy trigger based on a 1 GeV 
threshold for an ECL total energy sum [18,30]. The trigger efficiency depends 
on the average transverse momentum of the two tracks in the laboratory frame, 
= p'P = [{yV^^/2)‘^ — ■ (1 — I cos 6 **p)^/^, where the latter is the 

transverse momentum of p(p) in the 77 c.m. frame. We determine the trigger 
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efficiency as a function of pP, since each of the two-dimensional bins in 
and I cos 0 * 1 , where the number of events is measured, is associated with a pP 
value using the relation above. From the data, the trigger efficiency is 0.83 ± 
0.02 at pj^ = 0.55 GeV/c and 0.95 ± 0.05 at pj^ = 0.95 GeV/c. Gorrections 
for the Monte Garlo trigger efficiency are implemented according to the data, 
with a systematic error within 5%. 

The accuracy of the PID efficiencies has been checked by comparing Monte 
Garlo estimates to those based on data. The efficiency associated with each 
of the four PID conditions (Section 2) is studied, using events passing all 
selection criteria except the condition in question. The overall PID efficiency 
is ~ 92%, ~ 88 % and down to ~ 78% at = 2, 3 and 4 GeV respectively, 
with a systematic error less than 6 % in the whole range. The fake rate is 
~ 0.01% - 0.3% for = 3-4 GeV, respectively. 

Monte Garlo studies indicate that the PID requirements are very efficient in 
rejecting electrons and other relativistic particles, so that events from 77 — 
7r~^7r~, fj,~ and do not survive the selection, leaving those from 77 — 
K^K~ as the main residual background. From Monte Garlo simulation and 
the measured cross sections for 77 — K^K~ [27,31], the fraction of data 
that can be attributed to residual K~^K~ background, fm, is evaluated. Based 
on Monte Garlo studies, the dependence of fm on | cos 9* \ is negligible and 
= 0.8 ± 0.3%, 3.2 ± 0.5% and 7.7 ± 0.8% at = 3.2, 3.6 and 4.0 
GeV, respectively. For IF^.y < 3.0 GeV, fm is negligible. The values of fm have 
been checked in the data, using events passing all selection criteria except 
that on the normalized likelihood, Xp > 0.8. The number of signal events 
that would pass all selection criteria is estimated from the AT distribution 
of one of the two tracks, after requiring the other to satisfy Xp > 0.8. The 
values of fm inferred in this way are in good agreement with those above. The 
contribution from this source of background is subtracted from the data, using 
the expected fm from the Monte Garlo studies. The systematic uncertainty 
due to this source is ~ 1 % or less in the whole range. 

Possible non-exclusive backgrounds {ppX), most of them from 77 — ppir^ 
events, have been searched for in the data. Monte Garlo studies show that a 
high purity sample of such background can be derived from events with larger 
iSpJ'l and smaller |SpJ‘(pp7r°)|, the latter being the transverse momentum bal¬ 
ance of the three particles. By comparing Monte Garlo and data distributions 
of these parameters we obtain the fraction of the data attributed to this back¬ 
ground type, fn- We hnd that the dependence of /„ on | cos0*| is negligible, 
and it ranges from 5 ± 2 % to 17 ± 8 % for from 2 to 4 GeV, respectively. 
Gorrections are made using the fPW^P obtained above, and in total 7±1% of 
the selected data are subtracted. The systematic error from this source is 2 - 
12% for from 2 to 4 GeV, respectively. The fraction of the data attributed 
to 77 — ppK^ events is also obtained as a function of iSpJ'l. Before the cor- 


14 



rection, a comparison of the | | distribution between data and Monte Carlo 

exhibits a total difference of ~ 9% between the two samples, while it is re¬ 
duced to less than 3% for any range after the correction (Fig. 7). Since the 
residual excess in the iSpJ'l distribution could be attributed to residual non¬ 
exclusive backgrounds and a broader nature of the signal distribution than the 
Monte Carlo, the systematic uncertainty due to other possible non-exclusive 
backgrounds is limited to 3% after the correction. 


W^= 2.0-4.0 GeV = 2.0 - 4.0 GeV 



= 2.0 - 2.4 GeV = 2.0 - 2.4 GeV 




Fig. 7. \T,pl\ distributions for the data before (left column) and after (right column) 
the subtraction of residual non-exclusive backgrounds (77 —> ppir^). The Monte 
Carlo distributions are scaled with the first bin normalized to the data. 


6 Theoretical approaches 


From the asymptotic QCD prediction of Eq.(l) and after integration over 
cos^*, the cross section for 77 —> pp is proportional to for asymptot¬ 

ically large hFy.y. All models based on this framework behave asymptotically 
as cr oc . For the diquark scenario, two curves are provided [14]: from 
the complete diquark model and from the same model with only helicity con¬ 
served amplitudes, where p and p are in opposite helicity states. The scale of 
the diquark predictions matches the data for = 2.5 - 3.0 GeV, but the 
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Belle (IcosG I < 0.6) Belle (IcosO I < 0.6) 



(GeV) (GeV) 


Fig. 8. Separate fits of ct oc W~^ to the data in the range of = 2.5 - 2.9 GeV 
and 3.2 - 4.0 GeV, with (a) n floating; (b) n = 10 and n = 15. The error bars include 
statistical and systematic errors. The x^/ndf values for each fit are indicated in the 
figure. 

deviation becomes larger as increases (Fig. 4). At higher energies, the 
data fall below the diquark predictions and exhibit a gradual approach to the 
three-quark model predictions [4]. At medium energies between 2.5 and 4.0 
GeV, a steeper fall of the total cross section in is observed. 

If we £t the data with a power law a oc W~^ with n floating (Fig. 8 (a)), 
taking into account both statistical and systematic uncertainties as well as 
possible correlations between the latter, we obtain n = 15. it?;® and I 2 . 4 II 3 
in the range of = 2.5 - 2.9 GeV and 3.2 - 4.0 GeV, respectively (the 
charmonium region between 2.9 and 3.2 GeV is excluded). For completeness, 
we also show in Fig. 8 (b) the results of the hts with n hxed at 10 and 15. 
Although for both ranges a good fit to the data can be obtained at n = 15, 
a smaller power, n = 10, describes the data above 3.2 GeV reasonably well. 
This may imply that lower power terms become dominant at higher enegies, 
which is an indication for the transition to the asymptotic predictions. 

The angular differential cross section in | cos 9* \ is another observable most 
important to the study of the nature of the interactions involved in the 
process 77 — pp. All existing models based on the constituent scattering 
picture [4,5,12,13,14,15], as expected, predict an ascending trend, which is 
in agreement with the data for > 2.5 GeV. This is due to the fac¬ 
tor l/V^u cc 1 /Vl — cos^ 6* contained in the hard scattering amplitudes. 
The same trend is obtained from naive QED [32] estimates: da/d\cos9*\ oc 
(1 -f- cos^0*)/(l — cos^0*), in the massless limit. A simplified picture with 
diquarks would follow the naive QED expectation above [10,12], if all quark 
masses are neglected and only scalar diquarks are considered. In Fig. 9 we plot 
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the differential cross section normalized to that averaged within | cos0*| < 0.3, 
and compare various predictions to the data. We observe that the data rise 
more sharply in | cos0*| at higher energy (see also Fig. 5). In comparison, all 
current models predict a flatter trend in | cos 6* \. 


2.5 < W^< 3.0 GeV 

YY 



3 < < 4 GeV 

YY 



Fig. 9. Differential cross section as a function of | cos 0*|, normalized to that averaged 
within I cos 9*\ < 0.3, for the two higher ranges of The error bars are statistical 
only. Theoretical predictions are from [14] (diquark), [5] (three-quark) and [15] 
(handbag). 


The deviation of the leading term QCD calculations [4,5] from the data at 
= 2.5 - 4.0 GeV implies that power corrections are still signihcant at 
these intermediate energies. It is not surprising since the very threshold of pp 
production corresponds to hFy.y ~ 2 GeV. However, the diquark and hand¬ 
bag models [12,13,14,15] were developed in order to describe the intermediate 
energy region at the price of introducing model form factors, etc.. The dis¬ 
agreement of the data at hFy.y = 2.5 - 4.0 GeV with their predictions (see 
Fig. 4 and 9) obviously necessitates their improvement. 


The descending trend of the differential cross section in |cos6**| observed at 
low energies (Hly-y < 2.5 GeV) cannot be understood within the hard scatter¬ 
ing picture (Fig. 6(a)). In a recent study based on non-perturbative QGD sum 
rules [33], this trend was proposed as a general feature for hadron pair pro¬ 
duction from two-photon collisions. The behavior is very natural if low partial 
waves are involved. In Ref. [16] it was shown that even a simple model based 
on pole- and resonance-dynamics can reproduce this behavior. 
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7 Conclusion 


Using the Belle detector at the high-luminosity KEKB collider, the cross sec¬ 
tions for 77 —pp have been measured for from 2.025 to 4.0 GeV and 
I cos 0*1 < 0.6, with systematic uncertainties ranging from 7% to 14%. These 
results represent a great improvement in precision compared to all previous 
measurements and allow more accurate tests of various theoretical models. We 
also observed the production of 77 —?7c —^ PP and determined the product of 
the two-photon width of the r]c and its branching ratio to pp. 


Fitting to a power law a oc shows that the best £t value of n decreases as 
energy increases, and n = 10 cannot be rejected at energies above 3.2 GeV, im¬ 
plying the gradual transition to the expectation from asymptotic predictions. 
The ascending trend for the differential cross section in | cos 6* \ predicted by 
the hard scattering picture is in agreement with the data for Wy.y > 2.5 GeV; 
however, the data rise more sharply in | cos 6* \ as increases. The descend¬ 
ing trend in |cos0*| at lower energies 14%^ < 2.5 GeV can be reproduced 
by non-perturbative approaches [16]. The descending trend of the differential 
cross section in | cos 6* \ changes to an ascending one with the increase of en¬ 
ergy, which could be an indication for the transition from a soft resonance 
regime to the beginning of a hard regime. Existing models suggested for the 
intermediate energies [12,13,14,15] can not provide satisfactory description of 
the observed energy and angular dependence in the studied energy range. 
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